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An  Analytic/Empirical  Model  of  the  Middle 
and  Low  Latitude  Ionosphere 

1.  INTRODUCTION 

Recently,  a  Semi-empirical  Low-latitude  Ionospheric  Model  (SLIM)  was  developed1  that  provides 
F-region  electron  and  ion  (0+)  density  profiles  from  180  km  to  1000  km  every  2°  latitude  (24°N  to  24°S) 
and  eveiy  one-half  hour  local  time  over  the  24-hour  day.  The  reason  for  producing  this  model  was  to 
generate,  in  a  computationally  efficient  manner,  electron  density  profiles  at  low  latitudes  that  were 
more  realistic  than  those  of  the  widely  used  empirical  models  of  Chiu2  or  Bent.3  SLIM  consists  of 
functions  and  coefficients  that  reproduce  9  sets  of  theoretically-calculated  low-latitude  electron 
density  profiles  corresponding  to  the  three  seasons:  equinox,  June  solstice,  and  December  solstice  for 
three  levels  of  solar  activity:  solar  minimum  (F10  7  =  70),  moderate  solar  activity  (F10  7  =  125),  and 
solar  maximum  (F 10  7  =  185).  A  separate  set  of  six  coefficients  exist  for  each  2°  interval  in  latitude  and 
each  one-half  hour  interval  in  local  time.  The  profiles  are  determined  by  the  six  parameters,  nmF2, 
hmF2,  HUP,  Hjjq,  Cup.  and  Clq  and  the  modified  Chapman  expression: 

Ne(Z)  =  nmF2  x  exp[C(l  -  Z  -  exp(-Z)]  (1) 


(Received  for  Publication  30  March  1989) 

1.  Anderson,  D.N.,  Mendlllo,  M.,  and  Hemiter,  B.  (1987a)  A  semi-empirical,  low-latitude 

ionospht  i  ic  model.  Radio  Set  22:292-306. 

2.  Chiu,  Y.T.  (1975)  An  improved  phenomenological  model  of  ionospheric  density.  J.  Atmos.  Terr. 

37:1563-1570. 

3.  Llewellyn,  S.K.  and  Bent,  R.B.  (1973)  Documentation  and  Description  of  the  Bent  Ionospheric 

Model.  AFCRL-TR-73-0657,  AD  772733,  Air  Force  Cambridge  Research  Laboratories,  Bedford, 
MA. 
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where  nmF2  is  the  peak  plasma  density,  hmF2  is  the  altitude  where  Ne(Z)  =  nmF2,  Z  =  h~1^|nF2i  and  the 
subscripts  UP  and  uq  refer,  respectively,  to  those  regions  above  and  below  hmF2.  While  electron 
density  profiles  are  easily  generated  from  the  various  sets  of  coefficients  given  in  tabular  form,  a  more 
convenient  approach  would  be  to  specify  the  coefficients  as  analytic  functions  of  local  time  and 
latitude,  as  in  the  Chiu2  empirical  model.  Furthermore,  because  the  Chiu  ionospheric  model  has  been 
so  extensively  used  to  specify  the  low-  and  middle-latitude  ionospheric  region,  it  was  decided  to  recast 
the  coefficients  in  an  analytic  form  so  that  present  users  of  the  Chiu  model  would  have  access  to  the 
new,  improved  formulation  without  having  to  modify  any  of  the  input  and  output  parameters. 

This  report  describes  briefly  the  technique  used  to  generate  the  analytic  functions  that  would 
reproduce  the  SLIM  profiles,  and  then  compares  this  Fully  Analytic  Ionospheric  Model  (hereafter 
referred  to  as  FAIM)  with  SLIM  and  the  Chiu  empirical  ionospheric  model  to  illustrate  the 
similarities  and  differences,  respectively.  We  then  demonstrate  the  importance  of  using  a  realistic 
low-latitude  ionospheric  model  when  calculating  the  local  time  variation  in  zonal  winds  from  the 
dynamical  neutral  wind  equations. 


2.  THE  APPROACH 


The  empirical  Chiu  ionospheric  model  is  based  on  middle-  and  low-latitude  observations  of  the 
E,  FI,  and  F2  regions  of  the  ionosphere  which  cover  seasonal  as  well  as  solar  cycle  variations.  The 
model  assumes  that  the  F2-region  plasma  density  can  be  expressed  by  the  Chapman  profile  expression 

Ne(Z)  =  nmF2  x  exp[l  -  Z  -  exp(-Z)),  (2) 


where  the  variables  are  the  same  as  in  Eq.  (1),  and  the  coefficient  C  is  taken  to  be  unity.  As  in  SLIM, 
the  Chiu  formulation  utilizes  two  expressions  for  Z,  one  for  the  bottomside  of  the  F-layer  and  one  for 
the  topside,  given  respectively  by 


h  -  hmF2 

Zlo  " 


(3) 


„  h  -  hmF2 
ZUP  -  Hup 


(4) 


In  all,  four  coefficients  are  needed  to  completely  specify  the  F2-layer  in  the  Chiu  model:  nmF2,  hmF2, 
HUP  and  H^.  Similar  formulations  are  utilized  for  the  FI-  and  E-layers,  and  the  total  plasma  density 
as  a  function  of  altitude  is  obtained  by  summing  these  three  contributions.  Analytic  functions  are 
used  throughout  to  express  the  coefficients  as  a  function  of  geographic  and  geomagnetic  coordinates, 
sunspot  number,  local  time,  and  day  of  year. 

The  Chiu  model  is  widely  known  to  be  deficient  in  its  specification  of  hmF2  at  middle  and  low 
latitudes,  and  its  representation  of  the  equatorial  anomaly.  In  our  approach,  we  have  first  applied  a 
correction  to  the  Chiu  hmF2s  at  middle  and  low  latitudes  based  on  ionosonde  data.  Then,  we  have 
added  analytic  corrections  for  foF2  and  hmF2  and  modified  HUP  and  Hlq  in  the  Chiu  model  between 
±24°  dip  latitude,  such  that  the  ionization  profiles  best  fit  those  of  the  SLIM  model.  In  FAIM  we  no 
longer  use  the  modified  Chapman  function  (1)  which  contains  the  CUP  and  Clq  coefficients.  The  above 
procedure  and  the  equations  involved  are  detailed  below. 
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3.  ANALYTIC  FORMULATION 


3.1  Middle  Latitude  Correction 

The  middle  latitude  correction  to  the  hmF2  specification  of  the  Chiu  model  is  based  on 
examination  of  a  rather  limited  set  of  monthly  mean  ionosonde  measurements  (about  40  stations) 
covering  only  September  1984  (a  solar  minimum  period  with  F10  7  =  80)  and  March  1979  (a  solar 
maximum  period  with  F10  7  =  180).  These  months  were  chosen  because  they  include  the  Equinox 
Transition  Study  (ETS)  campaign  and  the  CDAW-6  interval,  and  hence  the  data  were  readily  available 
to  us.  Our  main  interest  here  is  to  provide  a  significantly  improved  and  easy  to  use  analytic/empirical 
model  of  the  low-latitude  ionosphere.  However,  when  we  realized  how  poor  the  F-layer  peak  heights 
(hmF2s)  in  the  Chiu2  model  really  were,  we  decided  to  apply  a  first-order  correction  to  these,  prior  to 
deriving  the  low-latitude  correction  based  on  SLIM.  Certainly,  these  middle-latitude  corrections 
ought  to  be  scrutinized  and  improved  upon  using  data  with  broader  coverage  in  season  and  solar  cycle, 
but  this  effort  is  viewed  out  of  the  scope  of  the  present  study.  We  did,  however,  find  the  Chiu  values  for 
foF2  to  match  the  ionosonde  data  reasonably  well,  and  did  not  attempt  any  corresponding  corrections 
to  this  ionospheric  parameter. 

The  original  expression  for  hmF2  in  Chiu2  is: 

hmF2  =  240  +  75p  +  83p£cosA.m  +  30cos((|>  -  4.5|  Xml  -  rt)  +  10coskmcos^  (t  -  4.5)  (5) 

where 

t  =  annual  time  (days  of  year)  In  units  of  months  as  measured  from  15  December  of  the 
previous  year 

p  =  R/100 

R  =  the  monthly  smooth  Zurich  sunspot  number 
=  geomagnetic  latitude  (rad) 

C  =  sin8slnXm  where  8  =  solar  angle  (rad) 

<t>  =  local  time  angle  (rad)  measured  from  local  midnight. 

Note  that  the  local  time  variation  consists  only  of  a  diurnal  (24- hour  harmonic)  that  is  independent  of 
season  and  solar  cycle.  Our  correction  concerns  the  second  (solar  cycle)  and  fourth  (local  time)  terms 
in  the  above  equation.  The  former  is  replaced  by 

100  Vp  -  10.  (6) 

The  local  time  term  is  replaced  by 

(56  -  12Vp)cos(<j)  +  Xa)  +  (9  +  12Vp)cos(2<J)  -  4)  +  12cos(3<(>  -  ^p)  (7) 

where 

o-  f  Xm  ^  2 

Xa=  -2tc  +-^-e  v349j  .  (8) 

Thus,  we  have  modified  the  mean  value  and  solar  cycle  dependence  of  the  diurnal  average 
component  of  hmF2,  have  modified  the  diurnal  (24-hour)  component  to  reflect  a  solar  cycle 
dependence,  and  in  addition,  have  added  semidiurnal  (12-hour)  and  terdiumal  (8-hour)  components  to 
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the  expression  for  hmF2.  These  dependences  were  derived  by  examining  the  local  time  Fourier 
components  of  hmF2  and  foF2  from  global  ionosonde  data  during  September,  1984,  and  March,  1979, 
as  mentioned  previously.  These  data  are  illustrated  in  the  Appendixes  of  this  report  in  Figures  Al- 
A14  and  B1  -  B14,  respectively,  and  are  compared  there  to  values  from  the  expressions  above,  as  well  as 
to  the  low-latitude  corrections  derived  from  the  SLIM  model,  to  be  described  below. 

Comparisons  between  the  local  time  variations  of  hmF2  as  observed  at  Wakkanai  (45.4° 
latitude),  Dakar  (14.8°  latitude),  and  Moscow  (55.5°  latitude)  are  illustrated  in  Figures  1  and  2.  These 
are  typical  of  comparisons  with  most  other  stations.  In  Figure  1 ,  the  solar  cycle  dependence  or  hmF2 
at  Wakkanai  is  seen  to  be  well  modelled  by  the  above  equations,  and  represents  a  considerably  better 
fit  to  the  data  than  Chiu.2  Figure  2  illustrates  the  latitude  dependence  of  the  local  time  variation  of 
hmF2  during  solar  minimum  conditions.  Again,  the  newly  defined  hmF2  representation  reproduces 
the  ionosonde  data  extremely  well,  reflecting  an  impressive  improvement  over  the  Chiu  formulation. 

3.2  Low  Latitude  Correction 

In  developing  the  low-latitude  poi ..  .n  of  FAIM,  we  derive  new  analytic  formulations  for  the 
parameters  foF2,  hmF2,  HUP,  and  Hlq  in  the  Chiu  model  based  on  least-squares  fits  to  the  SLIM 
profiles.  The  local  time  variation  at  a  given  dip  latitude  (Xm)  is  expressed  as  a  series  of  sines  and 
cosines  with  amplitudes  an  and  bn  with  the  series  going  from  n  =  0  to  n  =  6: 

/(X,<|>)  =  V6  (an(A.)sin(n<)>)  +  bn(X)cos(n<j>))  (9) 

■fc—n=0 

where  <J>  is  local  time  from  midnight  in  radians.  A  separate  function  fs(X,  4>)  is  found  for  each  of  the 
four  defining  parameters,  s  =  1  to  4. 

In  reality,  for  foF2  and  hmF2  values,  we  used  the  expression  above  to  fit  the  difference  between 
the  Chiu  and  SLIM  values  so  that  in  the  new  analytic  model,  we  simply  calculate  the  Chiu  values  and 
then  add  to  these  the  analytically-formulated  difference  values  to  provide  the  final  foF2  and  hmF2 
quantities  specified  by  FAIM.  The  procedure  for  obtaining  the  topside  and  bottomside  F-region  scale 
heights,  HUP  and  HIX)  is  handled  somewhat  differently.  Here,  a  separate  topside  and  bottomside  scale 
height  was  found  which  would  reproduce  each  SLIM  profile  density  to  within  0.5  MHz  assuming  the 
expression  for  Ne(Z)  given  by  Eq.  (2).  These  scale  heights  are  different  from  those  specified  in  the  paper 
by  Anderson  et  al1  because  the  values  for  Cyp  and  Cyo  are  now  assumed  to  be  unity.  Once  values  of  Hyp 
and  H[X)  were  obtained  by  this  method,  the  analytic  functions  given  by  Eq.  (9)  were  derived. 

The  next  task  was  to  find  appropriate  functions  and  coefficients  that  would  reproduce  the 
separate  sets  of  an's  and  bn’s  as  a  function  of  dip  latitude  (Xm).  The  modified  harmonic  oscillator 
function  Hk(y)e~'*'2  where  Hk  is  the  Hermite  Polynomial  of  order  k  and  y  =  ^pj  in  radians)  was 
found  to  be  ideally  suited  for  this  purpose: 

anUm).bn(*m)  =  Zk=0CkHk(v)e"v2-  (10) 

The  FAIM  FORTRAN  code  is  basically  the  Chiu2  code  with  the  corrections  outlined  above 
incorporated. 
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Figure  2.  Monthly  Mean  hmF2s  at  Dakar  (14.8°N.  342.6°E  Geographic;  21.4°N,  56.0°E  Geomagnetic)  and  Moscow  (55.5C 
39.3°E  Geographic;  50.4°N,  123°E  Geomagnetic)  Compared  with  Chiu  and  FAIM  Values.  The  hmF2  values  are  computed 
Figure  1. 


4.  RESULTS 


As  outlined  in  the  paper  by  Anderson  et  al1,  nine  separate  SLIM  cases  were  generated  covering 
three  levels  of  solar  activity:  low,  moderate,  and  high  for  each  of  three  seasons:  equinox,  June  solstice, 
and  December  solstice.  For  these  nine  cases,  sets  of  coefficients  to  the  analytic  functions  described 
above  were  determined  that  would  reproduce  the  local  time  and  latitude  variations  in  the  four 
parameters  foF2.  hmF2,  HUP,  and  H^. 

To  demonstrate  how  well  FAIM  reproduces  the  SLIM  values,  we  plot  in  Figure  3  the  latitude 
variation  in  hmF2  and  foF2  at  1400  LT  for  the  solar  cycle  maximum,  equinox  case.  In  this  figure,  we 
also  compare  these  models  with  the  empirical  Chiu  values  for  hmF2  and  foF2.  Looking  first  at  the 
hmF2  comparisons,  FAIM  smooths  out  the  SLIM  values  between  12°N  and  12°S  dip  latitude  and  gives 
excellent  agreement  poleward  of  ±12°.  Actually,  the  large  variation  in  SLIM  hmF2  between  ±4°  and 
±10°  dip  latitude  is  somewhat  an  artifact  of  the  SLIM  least-squares  approach  so  that  the  smoother 
variation  in  FAIM  hmF2  values  is  actually  more  realistic.  The  comparison  of  both  SLIM  and  FAIM 
with  the  Chiu  model  illustrates  the  magnitude  of  the  improvement  at  low  latitudes.  As  reported  by 
Rajaram  and  Rastogi5,  daytime  hmF2  values  near  the  magnetic  equator  at  Kodaikanal.  India  lie  in  the 
range  450  to  500  km,  about  100  km  greater  than  the  350  km  specified  in  the  Chiu  model. 

The  comparison  in  foF2  between  FAIM  and  SLIM  is  excellent.  The  analytic  model  slightly 
overestimates  SLIM  values  near  the  crests  of  the  equatorial  anomaly  at  ±16°  to  ±20°  dip  latitude  and 
underestimates  slightly  near  the  magnetic  equator,  but  these  differences  are  well  within  the  day-to- 
day  fluctuations  in  equatorial  foF2  values.  A  comparison  with  the  Chiu  model  shows  the  crest-to- 
trough  density  ratio  is  greater  for  the  analytic  model,  which  is  in  good  agreement  with  observed  values 
in  the  Indian  sector  during  March  1959,  a  solar  cycle  maximum  equinoctial  period.5 

At  2000  LT,  the  comparison  between  SLIM,  FAIM,  and  Chiu  are  even  more  Impressive  as  shown 
in  Figure  4.  The  post-sunset  enhancement  in  upward  E  x  B  drift  is  a  characteristic  feature  during 
equinox,  solar  maximum  years6  and  is  responsible  for  the  very  high  hmF2  values  near  the  magnetic 
equator.7  It  is  reproduced  in  both  SLIM  and  FAIM.  Such  high  hmF2  values  are  also  in  agreement  with 
measured  values  at  Kodaikanal4  and  at  Jicamarca. 1  The  analytic  model  reproduces  the  latitude 
variation  in  SLIM  values  very  well  and  both  these  models  give  substantially  greater  hmF2  values  than 
the  Chiu  model. 

An  unrealistically  low  F-layer  peak  height  is  provided  by  the  Chiu  model  and  can  have 
significant  consequences  if  it  is  used  in  neutral  dynamics  calculations.  In  particular,  after  sunset 
when  the  two  main  horizontal  forces  on  the  neutrals  are  the  pressure  gradiant  and  ion  drag  forces,  the 
low  hmF2's  and.  hence,  unrealistically  large  plasma  densities  are  expected  to  result  in  significant 


5.  Rajaram,  G.  and  Rastogi,  R.G.  (1977)  Equatorial  electron  densities-seasonal  and  solar  cycle 

changes.  J.  Atmos.  Terr.  Phys.  39: 1 175- 1182. 

6.  Rastogi,  RG.  (1966)  The  equatorial  anomaly  in  the  F2  region  of  the  ionosphere,  Insf.  Telecomm. 

Engrs.  12:245-256. 

7.  Woodman,  R.F.  (1970)  Verticle  drift  velocities  and  east-west  electric  fields  at  the  magnetic 

equator,  J.  Geophys.  Res.  75:6249-6259. 
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Figure  3.  Comparison  of  hmF2  and  foF2  Values  as  a  Function  of  Dip  Latitude  at  1 400  LT  Given  by 
SLIM,  Chiu,  and  FAIM  for  Solar  Cycle  Maximum,  Equinoctial  Conditions. 
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Figure  4.  Comparison  of  hmF2  and  foF2  Values  as  a  Function  of  Dip  Latitude  at  2000  LT  Given  by 
SLIM.  Chiu,  and  FAIM  for  Solar  Cycle  Maximum.  Equinoctial  Conditions. 
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underestimates  of  the  zonal  wind  speed.8 

The  enhancement  in  the  upward  E  x  B  drift  velocity  at  the  equator  mentioned  above  is  also 
responsible  for  maintaining  the  crests  of  the  anomaly  out  to  ±18°  dip  latitude  and  increasing  the  peak- 
to-trough  ratio  after  sunset.  This  latitude  variation  in  foF2  at  2000  LT  is  completely  absent  in  the 
Chiu  model,  which  essentially  displays  a  constant  density  value  with  latitude.  The  more  realistic 
FAIM  and  SLIM  models  with  peak  densities  of  5  x  106  el/cm3  (foF2  ~  18  MHz)  at  the  crests  of  the 
anomaly,  provide  much  greater  ion  drag  forces  at  these  latitudes,  which  translate  to  a  greater  latitude 
variation  in  the  zonal  neutral  wind  velocities.  In  addition,  this  latitude  structure  in  peak  densities 
significantly  affects  flux  tube  integrated  conductivities  after  sunset  that  play  an  important  role  in  the 
electrical  coupling  between  the  low-latitude  ionosphere  and  neutral  atmosphere.9 

A  comparison  of  profile  shapes  at  the  magnetic  equator  and  at  16’S  dip  latitude  is  illustrated  in 
Figure  5  for  the  2000  LT  period.  At  the  magnetic  equator,  reasonable  agreement  between  FAIM  and 
SLIM  is  achieved  with  the  topside  portion  being  in  better  agreement  than  the  bottomside  region. 
Actually,  this  represents  one  of  the  poorest  fits  of  all  the  comparisons.  The  degree  to  which  both  FAIM 
and  SLIM  profiles  differ  from  the  Chiu  model  is  graphically  illustrated  and  displays  the 
unrealistically  large  electron  density  values  in  the  350  km  region  provided  in  the  Chiu  model.  At  16°S 
dip  latitude,  the  comparison  between  FAIM  and  SLIM  is  excellent  and  is  more  indicative  of  the  general 
agreement  that  can  be  achieved  when  analytic  functions  for  Hyp  and  H1X)  are  used  to  calculate  profile 
shapes.  Here,  FAIM  and  SLIM  foF2  and  hmF2  values  are  significantly  greater  than  Chiu  values. 


5.  NEUTRAL  DYNAMICS  IMPLICATIONS 


As  shown  by  Anderson  and  Roble8,  the  spatial  and  temporal  variability  of  F-layer  densities  near 
the  magnetic  equator  can  have  dramatic  effects  on  the  neutral  dynamics.  We  seek  here  to  demonstrate 
typical  differences  in  neutral  wind  structures  implied  by  the  more  realistic  FAIM  than  the  widely  used 
Chiu2  empirical  ionospheric  model.  To  do  this  requires  only  a  fairly  simple  model.  At  the  geographic 
equator,  the  Coriolis  force  is  zero  and  the  following  zonal  momentum  equation  decouples  from  the 
meridional  equation: 


5u 
5t  = 


—  §£ 
p  5x 


-  vp(u  -  up)  + 


LL 

P 


82u 

8z2 


where 

u  =  zonal  (eastward)  wind  speed 
t  =  local  time 
p  =  neutral  mass  density 
p  =  neutral  pressure 


(ID 


8.  Anderson,  D.N.  and  Roble,  R.G.  (1974)  The  effect  of  verticle  E  x  B  ionospheric  drifts  on  F-region 

neutral  winds  in  the  low-latitude  thermosphere,  J.  Geophys.  Res.  79:5231-5236. 

9.  Anderson,  D.N.,  Heelis,  R.A.,  and  McClure,  J.P.  (1987b)  Calculated  nighttime  eastward  plasma 

drift  velocities  at  low  latitudes  and  their  solar  cycle  dependence,  Annul.  Geophys.  5A(6):435- 
442. 
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Figure  5.  (a)  Comparison  Between  SLIM,  Chiu,  and  FAIM  Electron  Density  Profiles  at  the  Magnetic 
Equator  at  2000  LT.  (b)  Same  as  (a).  Except  at  16°S  Dip  Latitude. 


x  =  distance  in  zonal  direction  =  flat 
Q  =  earth's  angular  velocity 
a  =  earth's  radius 
p  =  molecular  viscosity  coefficient 
z  =  altitude 

vp  =  plasma- neutral  collision  frequency 
up  =  plasma  drift  speed 

Pressure  gradients  are  calculated  from  the  MSIS-83  model10  which  is  based  heavily  on  AE-E  data  at 
low  latitudes,  and  so  should  represent  pressure  gradients  with  reasonable  accuracy,  at  least  at  low  to 
moderate  solar  activity.  A  simple  diurnal  variation  of  zonal  plasma  drifts  is  assumed  based  on  data 
from  Jicamarca  radar1 1 


up  =  30  +  70  cos (t  -  22)  (12) 

in  ms-1  where  t  =  local  time  in  hours. 

Contours  of  eastward  wind  speed  as  a  function  of  altitude  (150  -  300  km)  and  local  time  (0  -  24 
hours)  are  illustrated  in  Figures  6a  and  6b  during  solar  maximum  conditions  (F10  7  =  185)  conditions 
for  the  Chiu  and  FAIM  models,  respectively.  Figure  6c  illustrates  the  differences  between  these  two 
contour  plots  (FAIM  minus  Chiu),  and  represents  the  differences  in  neutral  wind  speed  due  to  the 
differences  in  the  two  model  plasma  density  distributions.  Figure  4d  similarly  illustrates  these 
difference  fields  for  solar  minimum  (F10  7  =  70)  conditions.  One  immediately  notes,  from  the 
difference  field  contours,  that  the  dynamical  consequences  of  utilizing  FAIM  vs.  CHIU  only  become 
apparent  between  about  1800-0200  LT  at  SSMAX,  and  1300-2300  LT  for  SSMIN.  In  addition,  the 
difference  fields  are  positive  at  SSMAX  and  negative  at  SSMIN.  There  are  several  competing 
mechanisms  that  interplay  to  produce  these  results.  First,  the  only  significant  differences  in  foF2 
between  FAIM  and  SUM  occur  between  1200  LT  and  (1900  LT,  2200  LT)  at  (SSMIN.  SSMAX).  where  the 
Chiu  foF2s  exceed  those  of  SLIM  by  2  -  3  MHz.  Furthermore,  during  the  daytime  the  neutral  dynamics 
are  considerably  less  sensitive  to  differences  in  hmF2,  since  the  layer  is  very  broad  and  the  ion  drag 
effect  is  distributed  over  a  wide  range  of  altitudes.  Finally,  the  effects  of  ion  drag  are  considerably  less 
at  the  magnetic  equator  than  at  higher  latitudes,  since  vp(u  -  up)  is  much  smaller  than  vpu.  especially 
at  night  when  u  and  up  are  interactively  coupled  via  the  F-layer  dynamo  mechanism. 

Despite  the  above  ameliorating  factors,  the  difference  fields  shown  in  Figures  6c  and  6d  still 
represent  a  significant  fraction  of  the  total  zonal  wind.  During  both  SSMIN  and  SSMAX.  the  Chiu 
model  foF2s  exceed  those  of  FAIM  by  about  2-3  MHz  between  1200  and  2100  LT.  In  addition,  during 
SSMAX.  hmF2s  from  FAIM  so  far  exceed  (by  100  -  150  km)  those  from  Chiu2  that  the  F-layer  is 
essentially  lifted  up  out  of  the  region  below  300  km;  and  the  ion  drag  force  is  essentially  eliminated, 
allowing  u  to  accelerate  to  values  of  160  -  200  ms-1.  During  SSMIN  the  situation  is  quite  different.  The 
hmF2s  given  by  Chiu  and  FAIM  for  the  pre- /post-sunset  period  are  very  similar,  and  quite  low 


10.  Hedin,  A.E.  (1983)  A  revised  thermospheric  model  based  on  mass  spectrometer  and  incoherent 

scatter  data:  MSIS-83,  J.  Geophys.  Res.  88:10170-10188. 

11.  Fejer,  B.G.,  Kudeki,  W.,  and  Farley,  D.T.  (1981)  Equatorial  F-region  zonal  plasma  drifts, 

J.  Geophys.  Res.  90:12249-12225. 
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(-300  km)  as  compared  to  SSMAX  (~450  km).  And.  even  though  the  Chiu  foF2s  exceed  those  of  FAIM  by 
2  MHz.  the  much  broader  profile  shape  and.  hence,  larger  height-integrated  ion  drag  force  given  by 
FAIM  leads  to  the  negative  difference  fields  in  Figure  6d.  Similarly  important  dynamical  effects  occur 
near  the  anomaly  peaks  (roughly  ±15°  geomagnetic  latitude).  One  would  therefore  anticipate 
significant  differences  in  the  low-latitude  thermosphere  circulation  pattern  simulated  with  a  3-D 
general  circulation  model  when  using  the  more  realistic  FAIM  model  as  opposed  to  Chiu2,  or  in 
electrodynamic  simulations  of  E-/F-region  coupling  involving  field-aligned  integrated 
conductivities. 
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Appendix  A 


Latitude  Variation  in  the  Diurnal,  Semidiurnal  and  Terdiurnal  Phase  and 
Amplitude  Components  of  hmF2  and  foF2  in  Three  Longitude  Sectors 

during  September  1984 
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Table  1.  Ionosonde  Stations  Utilized  for  Data  Depicted  In  Figures  A1-A14 

September,  1984 


Station 

Kergulen 

Geomagne 

-57.4, 

:tie  Lat.  Lone 
+129.9 

La  Reunion 

-21.2,  +55.6 

-27.4, 

+121.5 

Tahiti 

-17.7,  +210.7 

-15.2, 

+284.4 

Manila 

+14.7,  +121.1 

+3.6, 

+  191.1 

Chung-Li 

+24.9,  +121.2 

+13.6, 

+  189.5 

Okinawa 

+26.3,  +127.8 

+15.5, 

-163.1 

Ouagadougou 

+12.4,  +358.5 

+16.2, 

+71.6 

Yamagawa 

+31.2,  +130.6 

+20.6, 

-160.9 

Maul 

+20.8,  +203.5 

+21.2, 

+269.6 

Dakar 

+14.8,  +342.6 

+21.4, 

+56.0 

Kokubungi-Tokyo 

+35.7,  +139.5 

+25.7, 

-153.3 

Akita 

+39.7,  +140.1 

+29.8, 

-153.2 

Tashkent 

+41.3,  +69.6 

+32.3, 

+145.2 

Alma  Ata 

+43.3,  +76.9 

+33.5, 

+151.9 

Wakkanal 

+45.4,  +141.7 

+35.5, 

-152.7 

Khabarousk 

+48.5,  +135.1 

+38.1. 

-158.7 

Karaganda 

+49.8,  +73.1 

+40.3, 

+149.8 

Irkutsk 

+52.5,  +104.0 

+41.2, 

+175.5 

Rome 

+41.8.  +12.5 

+42.3, 

+93.2 

Pt.  Arguello 

+35.6,  +239.4 

+42.3. 

+302.4 

Novosibirsk 

+54.6,  +83.2 

+44.2. 

+158.9 

Bekescsaba 

+46.7,  +21.2 

+45.2, 

+103.2 

Tomsk 

+56.5,  +84.9 

+46.0, 

+160.6 

Kiev 

+50.5.  +30.5 

+47.1, 

+113.3 

Boulder 

+40.0.  +254.7 

+48.9. 

+318.2 

Poitiers 

+46.6,  +0.3 

+49.2, 

+83.0 

Wallops 

+37.8,  +284.5 

+49.2, 

+353.9 

Gorky 

+56.2,  +44.3 

+50.2. 

+127.7 

Moscow 

+55.5,  +39.3 

+50.4, 

+123.2 

Dourbes 

+50.1,  +4.6 

+51.7, 

+88.9 

Lannlon 

+48.8,  +356.5 

+52.0, 

+80.1 

Kaliningrad 

+54.7,  +20.6 

+53.0, 

+106.4 

Slough 

+51.5,  +359.4 

+54.0, 

+70.1 

Juliusruh 

+54.6,  +13.4 

+54.3, 

+99.7 

Leningrad 

+60.0,  +30.7 

+56.1, 

+118.3 

Ottawa 

+45.1,  +283.9 

+56.3, 

+352.7 

Salekhard 

+66.5,  +66.5 

+57.4, 

+149.7 

Nurmijarvi 

+60.5,  +24.6 

+57.7, 

+113.5 

Uppsala 

+59.8,  +17.6 

+58.3. 

+106.9 

Arkhangelsk 

+64.4,  +40.5 

+58.7, 

+129.1 

Lycksele 

+64.6,  +18.7 

+62.5, 

+111.7 

Loparskaya 

+68.0,  +33.0 

+63.1, 

+126.3 

Sodankyla 

+67.4,  +26.6 

+63.6, 

+120.8 

Klruna 

+67.8,  +20.4 

+65.1, 

+  116.4 

Godhavn 

+69.3,  +306.5 

+79.6, 

+34.7 

Resolute  Bay 

+74.7,  +83.2 

+83.2, 

+292.9 
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Figure  Al.  Diurnal  and  monthly  mean  hmF2's  from  ionosonde  stations  in  the  three  longitude  sectors. 
Sector  1:  Europe,  Africa.  W.  Asia;  Sector  2:  E.  Asia.  Pacific;  Sector  3:  America.  Greenland,  during 
September,  1984,  and  compared  to  the  SLIM  and  FAIM  models. 
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Figure  A3.  Same  as  Al,  Except  for  Semidiurnal  Amplitude. 
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Figure  A4.  Same  as  A1 ,  Except  for  Terdiumal  Amplitude. 


90. 


22 


PHASE  (Hour) 


hmF2 


Figure  A5.  Same  as  Al,  Except  for  Diurnal  Phase. 
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Figure  A6.  Same  as  Al.  Except  for  Semidiurnal  Phase. 
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Figure  A7.  Same  as  A1 ,  Except  for  Terdiuma)  Phase. 
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Figure  A9.  Same  as  A2.  Except  for  foF2. 
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Figure  A10.  Same  as  A3,  Except  for  foF2. 
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Figure  All.  Same  as  A4 .  Except  for  foF2 . 
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Figure  A 12.  Same  as  A5,  Except  for  foF2. 
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Figure  A13.  Same  as  A6,  Except  for  foF2. 
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Figure  A14.  Same  as  A7,  Except  for  foF2. 
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Appendix  B 


Latitude  Variation  in  the  Diurnal,  Semidiurnal  and  Terdiurnal  Phase  and 
Amplitude  Components  of  hmF2  and  foF2  in  Three  Longitude  Sectors 

during  March  1979 
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Table  2.  Ionosonde  Stations  Utilized  for  Data  Depicted  in  Figures  B1-B14 

March,  1979 


Station 

Geographic  Lat.  Long 

Bankok 

+13.8,  +100.6 

+2.6, +169.9 

Manila 

+14.7,  +121.1 

+3.6, +191.1 

Djibouti 

+  11.5,  +42.8 

+6.9, +114.9 

Chung-Li 

+24.9,  +121.2 

+13.6, +  189.5 

Okinawa 

+26.3,  +127.8 

+15.5,-163.1 

Ouagadougou 

+12.4,  +358.5 

+16.2,  +71.6 

Yamagawa 

+31.2,  +130.6 

+20.6,-160.9 

Maui 

+20.8,  +203.5 

+21.2, +269.6 

Kokubungi-Tokyo 

+35.7,  +139.5 

+25.7,-153.3 

Seoul 

+37.2,  +126.6 

+26.3, +195.0 

Mexico 

+19.3,  +260.6 

+29.0,  +328.4 

Akita 

+39.7,  +140.1 

+29.8,-153.2 

Wakkanai 

+45.4,  +141.7 

+35.5,-152.7 

White  Sands 

+32.3,  +253.5 

+41.2.  +318.6 

Rome 

+41.8,  +12.5 

+42.3,  +93.2 

Pt.  Arguello 

+35.6.  +239.4 

+42.3.  +302.4 

Tomsk 

+56.5,  +84.9 

+46.0. +160.6 

Boulder 

+40.0,  +254.7 

+48.9.  +318.2 

Poitiers 

+46.6,  +0.3 

+49.2,  +83.0 

Wallops 

+37.8,  +284.5 

+49.2,  +353.9 

Miedzeszyn 

+52.2,  +21.2 

+50.5. +105.7 

Yakutsk 

+62.0,  +129.6 

+51.2. +  194.9 

Dourbes 

+50.1,  +4.6 

+51.7,  +88.9 

Lindau 

+51.7,  +10.1 

+52.1,  +95.0 

Lannion 

+48.8,  +356.5 

+52.0.  +80.1 

Slough 

+51.5,  +359.4 

+54.0.  +70.1 

Juliusruh 

+54.6,  +13.4 

+54.3,  +99.7 

De  Bilt 

+52.1,  +5.2 

+53.5,  +90.5 

Ottawa 

+45.1,  +283.9 

+56.3.  +352.7 

Nurmljarvl 

+60.5,  +24.6 

+57.7, +113.5 

Uppsala 

+59.8,  +17.6 

+58.3. +106.9 

Lycksele 

+64.6,  +18.7 

+62.5. +11 1.7 

Sodankyla 

+67.4,  +26.6 

+63.6, +  120.8 

Kiruna 

+67.8,  +20.4 

+65.1, +116.4 

Resolute  Bay 

+74.7,  +83.2 

+83.2, +292.9 
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Figure  Bl.  Diurnal  and  monthly  mean  hmF2's  from  ionosonde  stations  in  the  three  longitude  sectors. 
Sector  1:  Europe,  Africa.  W.  Asia;  Sector  2:  E.  Asia,  Pacific;  Sector  3:  America.  Greenland,  during 
March.  1979,  and  compared  to  the  SLIM  and  FAIM  models. 
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Figure  B3.  Same  as  Bl,  Except  for  Semidiurnal  Amplitude. 
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Figure  B4.  Same  as  Bl,  Except  for  Terdiumal  Amplitude. 
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Figure  B5.  Same  as  Bl,  Except  for  Diurnal  Phase. 
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Figure  B6.  Same  as  Bl.  Except  for  Semidiurnal  Phase. 
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Figure  B7.  Same  as  B1 ,  Except  for  Terdiumal  Phase. 
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Figure  B8.  Same  as  B1 ,  Except  for  foF2. 
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Figure  B9.  Same  as  B2,  Except  for  foF2. 
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Figure  BIO.  Same  as  B3,  Except  for  foF2. 
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Figure  B1 1.  Same  as  B4,  Except  for  foF2. 
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Figure  B 12.  Same  as  B5,  Except  for  foF2. 
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Figure  B13.  Same  as  B6,  Except  for  foF2. 
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Figure  B14.  Same  as  B7,  Except  for  foF2. 
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